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We construct a model for the formation of black holes within galactic bulges. The 
\£> • initial state is a slowly rotating isothermal sphere, characterized by effective transport 

speed a e fr and rotation rate Q. The black hole mass is determined when the centrifugal 
radius of the collapse flow exceeds the capture radius of the central black hole. This 
model reproduces the observed correlation between black hole masses and galactic ve- 
t^J- \ locity dispersions, Mbh ~ 1O 8 M0(ct/2OO km s -1 ) 4 , where a = \/2a e g. This model also 

predicts the ratio \ib of black hole mass to host mass: hb ~ 0.004 (c/200 km s _1 ). 
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Two groups have recently reported an observed relationship between a galaxy's velocity disper- 
sion a and the mass Mbh of its central (supermassive) black hole. This correlation can be written 
in the form 

M bh = M (cr/200 km s" 1 ) 7 , (1) 



where the two groups find 7 = 3.75 (Gebhardt et al. 2000) and 7 = 4.7 (Farrarese & Merritt 
2000) and where the leading mass coefficient is Mo ~ 1.2 — 1.3 x 10 s Mq. Both groups report 
little scatter about the relationship (the scatter in Mbh at fixed a is bounded by 0.30 dex) and find 
no evidence for a subsidiary dependence on the Hubble type, profile type, or environment. This 
relationship supersedes a previous one claimed between the black hole mass and bulge luminosity 
(Richstone et al. 1998; Magorrian et al. 1998; van der Marel 2000; Kormendy & Richstone 1998). 
Furthermore, this measured scaling law poses a clear challenge for galaxy formation theories, which 
must ultimately account for this relationship. 

Several theories imply a relationship between the central black hole mass and the galactic 
velocity dispersion. For example, a semi-analytic model of merger-driven starbursts with black 
hole accretion (Haehnelt & Kauffman 2000; Kauffman & Haehnelt 2000) provides a correlation of 
the form [1]. Several models are based on black hole accretion influencing star formation and gas 
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dynamics in the host galaxy; this feedback can occur through ionization, mechanical work, and 
heating (e.g., Ciotti & Ostriker 1997, 2000; Blandford 1999; Silk & Rees 1998). The model of 
Blandford (1999) gives M bh < 77 a 5 / 3 whereas Silk & Rees (1998) predict M bh oc a . Finally, the 
accretion of collisional dark matter indicates the scaling relation Mbh cx a 4-4 ' 5 (Ostriker 2000). 

In this letter, we present a new theory of the Mbh — o relation, based on an idealized model 
for the collapse of the inner part of protogalaxies. The simplest variant of this theory is described 
in terms of particle dynamics in §2. The role of gas dynamics is explored in §3. We discuss the 
implications of this model and identify future issues in §4. 



2. THE BASIC BALLISTIC MODEL 

In this section, we examine the collapse of the inner part of a region that will form the bulge 
of a galaxy. The calculation starts at the time of maximum expansion for the main body of the 
bulge. For the sake of definiteness, we assume the following: [1] The dark matter and baryons 
are unsegregated. [2] The mass in this region is distributed like a singular isothermal sphere even 
though it is not in virial equilibrium. The initial density and mass distributions thus take the form 

n 2 In 2 
P{r) = 2^ and M(r) = -fr. (2) 

The transport speed a e g that specifies the initial conditions is related to the isotropic velocity 
dispersion a according to a = v^Oefr (see below). [3] This region is slowly rotating like a solid body 
(e.g., due to tidal torques) at a well-defined frequency Q. Both dark matter particles and parcels 
of baryons that are initially located at radius Too thus have initial angular momentum j = r^O 
sin 2 9, where 9 is the polar angle in spherical coordinates. [4] The center of this region contains a 
black hole, which may have a tiny mass at the start. This initial "seed" black hole could form by 
the collapse of the densest (central) part of the perturbation or could be primordial. 

Particles in the main body of the initial distribution will fall toward the center. Because the 
dynamical time scales monotonically increase with radius, infalling shells do not cross. The mass 
contained inside a given spherical shell, which marks a particle's location, does not change as the 
particle falls inward. As a result, orbital energy is conserved and is given (classically) by 

„ 1 2 1 f GM 
2 lr z r 

We consider orbits which fall a long way toward the center of the galaxy. Idealizing these trajectories 
as zero energy orbits permits the use of equation [3] to determine their pericenters p. For particles 
falling within the equatorial plane (9 = 7r/2), the pericenter can be written in the form 

P 2GM 2GM 2 5 a 8 eS ' 1 ' 



-3- 



where is the starting radius of an infalling particle and where we have used M = M(r) as a 
label for in the final equality. 

If this pericenter p is sufficiently small, ballistic particles will pass inside the horizon of the black 
hole and be captured. As mass accumulates in the black hole, its horizon scale grows accordingly. 
The pericenter of particles in ballistic orbits, falling from our assumed mass distribution, increases 
as p oc ~ M 3 . In the earliest stages of the collapse, all of the falling material is thus captured 
by the black hole. Later, this growth mechanism cuts off sharply when the black hole mass reaches 
a critical point defined by equating the pericenter p (for 6 = tt/2 orbits) to the capture radius of 
the black hole. In Schwarzschild geometry, particles coming in from infinity on zero energy orbits 
are captured by the black hole if p < AR$ (Misner, Thorne, & Wheeler 1973), where Rs = 2GM/c 2 
is the Schwarzschild radius. The condition p = ARs thus defines the critical mass scale Mc where 
direct accretion is compromised, i.e., 

In this scenario, the critical mass scale Mq determines the observed black hole mass Mbh (note 
that equation [5] displays the correct scaling with a). Most of the baryonic material not captured 
by the black hole during this early collapse phase eventually forms stars in the galactic bulge. Dark 
matter with low angular momentum is captured into the black hole along with the baryons; dark 
matter with high angular momentum (p > 4Rs) passes right through the galactic plane and forms 
an extended structure. 

We now use observations to specify the appropriate values of Q and a e g. In dissipationless 
collapse, the scale length of the mass distribution drops by 1/2 from maximum expansion, and the 
"flat rotation curve conspiracy" suggests that dissipation doesn't alter a further. This argument 
implies that the observed velocity dispersion a is related to the initial transport speed a c g of the 
protogalactic material through the relation a 2 = 2a 2 s . To specify f2, we use the fundamental 
plane, which provides a well defined relationship (for normal ellipticals and bulges) between the 
half-light radii of galactic bulges and the corresponding velocity dispersions (see Binney & Merrifield 
1998). For a = 200 km s _1 , the effective radius R of a mean surface brightness elliptical on the 
fundamental plane is about 3.5 kpc. With this scale R as the outer boundary of the forming bulge, 
the angular momentum of a ballistic particle cannot exceed the angular momentum of a circular 
orbit at this radius, i.e., R 2 Q. = oR or f2 = a/R. Evaluating this result for a = 200 km s _1 , we 
obtain a fiducial value of Q = 5.8 x 10~ 2 Myr _1 = 1.8 x 10~ 15 rad s _1 . With the values of Q and a 
now specified, we use equation [5] to find the desired Mbh — c relation 

M bh ~ 1O 8 M (cj/2OO km s" 1 ) 4 , (6) 



where we have written the result in terms of a rather than a e s . This relation is in good agreement 
with the observed correlations (equation [1]), as shown in Figure 1. 
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3. GAS DYNAMICS AND OTHER EFFECTS 

In this section, we include gas dynamics in our model for black hole formation during the 
collapse of galactic bulges. Because these structures become gravitationally bound, we expect the 
proto-bulge to collapse as a whole from an initial state (described here by equation [2]). To obtain 
a mathematical description of this collapse, we consider the flow that produces a galactic bulge to 
be a scaled up version of the collapse flows that have been studied previously for star formation 
(Shu 1977; Terebey, Shu, & Cassen 1984); this approach should thus capture the basic essence 
of the collapse problem. The collapse of the initial state (with density distribution [2]) proceeds 
from inside-out and the central portion of the flow approaches a ballistic (pressure-free) form: Dark 
matter always exhibits pressure-free behavior. Even for infalling gas, however, the inner limit of the 
collapse flow approaches pressure-free conditions. For collapse over sufficiently long time intervals, 
stars can form as parcels of gas fall inward and the resulting infalling stars are manifestly ballistic. 
The time scale for individual star formation events is r* ~ 10 5 years (e.g., Adams & Fatuzzo 1996), 
whereas the time scale Tbi g for the entire bulge structure to form is much longer (rbi g ~ 25 - 50 
Myr) . We thus expect most of the stars to form while the overall collapse of the bulge is still taking 
place. 

For a given gravitational potential, we find the orbital solutions for stars (or gas parcels or 
dark matter) falling towards the galactic center. In this initial calculation, the inner solution is 
derived using the gravitational potential of a point source. This form is only used in the innermost 
regime of the collapse flow where the potential is dominated by the forming black hole. As a result, 
this orbital solution is valid for the range of length scales Rs < r <C r^. 1 Since this potential 
is spherically symmetric, angular momentum is conserved and the motion is confined to a plane 
described by the coordinates (r, (f>); the radius r is the same in both the plane and the original 
spherical coordinates. The angular coordinate 4> in the plane is related to the angle in spherical 
coordinates by the relation cos (j) = cos / cos 6>o, where 6q is the angle of the asymptotically radial 
streamline (see below). For zero energy orbits, the equations of motion imply a cubic orbit solution, 

i - if. = (i - jig) -iL , (7) 

/io GMr 

where joo is the specific angular momentum of particles currently arriving at the galactic center 
along the equatorial plane. Here, the trajectory that is currently passing through the position (r, 
fi = cos 9) initially made the angle 6q with respect to the rotation axis (where = cos#o)- As in 
star formation theory, we define the centrifugal radius Rc = j^/GM (= 2p), which represents the 
radius of a circular orbit with angular momentum joo (for 6 = it/ 2). Our assumption of uniform 
initial rotation at rate f2 implies that joo = fJr^, where is the starting radius of the material 
that is arriving at the center at a given time. To evaluate the radii Rc and r^, we invert the mass 
distribution of the initial state (equation [2]) to find r^GM /2a 2 cS and R c = 2 G 3 M 3 /16a^ ff . 



1 Notice that at late times, long after black hole formation is complete, the potential is no longer close to a point 
potential and this solution loses its validity; in addition, relativistic corrections become important as r — ► Rs- 
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With an isothermal profile as the initial state, the collapse solution indicates that the flow- 
exhibits a well defined mass infall rate M = moa^/G, where rriQ ~ 0.975 (Shu 1977). This infall 
rate is constant in time and we can measure the time since the collapse began by the total mass 
M that has fallen to the galactic center. At early times, all of the mass falling to the center is 
incorporated into the central black hole. At later times, the mass supply is abruptly shut off by 
conservation of angular momentum. In this setting, the mass infall rate is quite large, M « 650 
M & yr -1 (for a = 200 km s _1 and a 2 = 2a^ ff ). The time scale T D h to form a typical supermassive 
black hole (with mass Mbh ~ 10 8 Mq) is thus about Tbh ~ 10 5 yr, comparable to the time scale 
for individual stars to form. The time scale to form the entire bulge is much longer, about 
Tbig ~ 25 — 50 Myr, comparable to the crossing time t croS s = R/o-cS- 

In this collapse flow, streamlines entering the central region do not cross each other. As long 
as the (bulge) infall time is longer than the time scale for individual star formation events, the core 
regions that produce stars will not interact. However, this scenario has an initial transient phase 
(~ 10 5 yr) in which the collapse of the bulge structure takes place faster than individual stars form. 
If stars are already condensing out of the collapse flow during this initial phase, they can interact 
and merge upon entering the central region. The resulting merger activity can lead to the initial 
production of the central black hole, which then gains additional mass as the collapse proceeds. 

Given the orbital solution (equation [7]), we can find the velocity fields for the collapse flow. 
The density distribution p(r, 6) of the infalling material can then be obtained by applying con- 
servation of mass along a streamtube (Terebey, Shu, & Cassen 1984) and can be written in the 
form 

( a\ _M_dno 
P{r > 6) = 4^pl^- (8) 

The properties of the collapsing structure determine the orbit equation [7], which in turn determines 
the form of d^o/dfj, and the radial velocity v r . The density field is thus completely specified 
(analytically, but implicitly). 

This gas dynamical version of the model defines the same critical mass scale for central black 
holes as the ballistic model presented in §2. In the earliest stages of collapse, incoming material 
falls to small radii r <C Rs, where Rs = 2GM/c 2 and M is the total mass M = Mt that has fallen 
thus far. In this early stage, the black hole mass Mbh = M. As the collapse develops, incoming 
material originates from ever larger radii and carries a commensurate increase in specific angular 
momentum. The centrifugal barrier of the collapse flow grows with time. The black hole mass is 
determined when the centrifugal radius exceeds the capture radius of a black hole in Schwarzschild 
geometry. This condition takes the form Rc > aR$, where a = 8 and is determined by particle 
orbits in a Schwarzschild metric. This condition, Rc > aRs, leads to the same critical mass scale 
Mq defined in equation [5]. Although Mc determines the black hole mass Mbh, two additional 
effects conspire to make the final black hole mass somewhat larger: 

[1] Even after the centrifugal barrier grows larger than the capture radius, the black hole 
continues to gain mass from infalling streamlines oriented along the rotational poles of the system. 
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The fraction of the infalling material that lands at such small radii is a rapidly decreasing function 
of time. As a result, this effect makes the black hole mass larger by only a modest factor Ta- The 
mass infall rate Mbh for material falling directly onto the black hole itself is given by 



M bh = j ^2ir(aRs) 2 \v r \p(aR s ,n), (9) 

where [i = cos 9. Using equation [8] to specify the density, we evaluate the integral to obtain a 
differential equation for the time evolution of the black hole mass. Solving the resulting differential 
equation (Adams et al. 2001), we find that the black hole mass increases by a factor Ta ~ 1-35 
due to direct infall. 

[2] Gas that falls to the midplane of the system can collect into a disk structure surrounding the 
nascent black hole. The presence of the disk is consistent with the current theoretical ideas about 
AGNs and the jets they produce. In order to retain the desired scaling law Mbh ~ c 4 , however, the 
total mass added to the black hole through disk accretion must be less than (or comparable to) the 
original mass scale Mq- We can estimate the maximum amount of mass that can be added to the 
black hole through the disk by using the constraint that disk accretion cannot operate faster than 
the orbit time at the outer disk edge (which is determined by Rc)- For reasonable assumptions, 
this maximum mass scale in the limit of efficient disk accretion is about 10 Mc (Adams et al. 
2001). 

This model also predicts a mass scale Mb for the bulge itself. If the initial protobulge structure 
is rotating at angular velocity Q, then only material within a length scale R = a e s/£l can collapse 
to form the bulge. Material at larger radii, r > R, is already rotationally supported and will not 
fall inwards. The length scale R thus defines an effective outer boundary to the collapsing region 
that forms the bulge. This boundary R, in turn, defines a mass scale for the bulge, i.e., 

2a 3 

M B = — j£ » 2.4 x 10 10 M Q (a/200 km s" 1 ) 3 . (10) 

Cxi I 

Both the bulge mass Mb and the black hole mass Mbh ~ Mc have the same dependence on the 
rotation rate We can divide out the rotation rate and find a robust estimate for the mass fraction 
[IB, i-e., 

H B = -J*L = ^32° ~ 0.0038 (cr/200 km s" 1 ). (11) 
Mb c 

This mass fraction jib is roughly comparable to the observed ratio of black hole masses to bulge 
masses in host galaxies (e.g., Richstone et al. 1998; Magorrian et al. 1998), although the data show 
appreciable scatter (see Figure 2). 



4. CONCLUSION 



In this paper, we have presented a simple model to describe the collapse flow that produces 
galactic bulges and the supermassive black holes living at their centers. The initial (pre-collapse) 
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state is a slowly rotating isothermal sphere, characterized by an effective sound speed a c g- and 
an angular velocity fi. These parameters (a e g,Q) thus represent the specification of the initial 
conditions. The velocity dispersion of the final stellar system is comparable to the initial sound 
speed and we make the identification a m v2a c ff- In developing this basic picture, we find the 
following results: 

[1] The black hole mass Mbh is determined by the condition that the centrifugal radius ex- 
ceeds the capture radius of a Schwarzschild black hole. This requirement leads to the scaling law 
Mbh = Mo(<r/200 km s -1 ) 4 , which is consistent with observations both in its dependence on veloc- 
ity dispersion a and for the mass scale (Mo ~ 10 8 M Q ) of the leading coefficient (see equation [6] 
and Figure 1). 

[2] The bulge mass is determined by the outer boundary of the collapsing region - material 
at initial radii r > R is rotationally supported and cannot collapse. This condition leads to the 
scaling law Mb oc cr 3 (equation [10]) and predicts the ratio //^ of black hole mass to bulge mass 
(equation [11]). This mass fraction scales weakly with velocity dispersion and has a typical value 
Hb ~ 0.004, similar to observed mass ratios (Figure 2). 

Elliptical galaxies (and bulges) are described by four parameters: effective radius R e , luminos- 
ity L, velocity dispersion a, and central black hole mass Mbh- For this discussion, we substitute 
the total mass Mb for the luminosity (for a given dark matter fraction, we thus assume that the 
transformation Mb —>■ L is determined by known, but perhaps complicated, stellar physics). We 
need to make the connection between these four basic quantities (a, R e , Mb, Mbh) and our model, 
which is described by only two variables: the transport speed a e g and the rotation rate f2. Our 
model implies the following transformation between initial conditions and final system properties: 
[1] a = V2a cS , [2] R e = a cS /(2fl), [3] Mb = 2a 3 eS /(Gn), [4] M bh = 16F A a 4 cS / \cGSl) , and [5] 
fiB = 8-7^4 Geff/cj where Ta ~ 1-35. The fifth relation is not independent, but because the rotation 
rate cancels out, it represents a robust prediction. The implications of these five scaling laws must 
be tested further against observations. 

In this letter, we have presented a working collapse model for black hole and bulge formation. 
In future work, a number of issues should be addressed, including relativistic corrections to the 
infall solutions, additional black hole mass contributions from disk accretion, and feedback from the 
central black hole on the collapse. In order to successfully reproduce the observed correlations, our 
model uses initial conditions with a particular profile of specific angular momentum (determined 
by the initial p(r) and Jl); we thus need to explore more general initial conditions and make 
the connection between the required initial conditions and earlier stages of galaxy formation. In 
addition, galactic bulges are likely to experience merger events (White & Rees 1978) which lead to 
modest increases in the overall velocity dispersion (White 1979); we must determine how the black 
holes produced through this model evolve in the context of such galactic mergers. In any case, the 
properties of these supermassive black holes - and their connection to their host galaxies - will 
ultimately provide a vital diagnostic for the galaxy formation process. 
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FIGURE CAPTIONS 

Figure 1. The correlation between black hole mass Mbh and velocity dispersion a of the host 
galaxy. The data points (adapted from Gebhardt et al. 2000) represent the observed correlation 
for ellipticals (circles), SO galaxies (squares), and spirals (triangles). The solid curve shows the 
theoretical result of this paper (using equation [6] and the 35% correction from equation [9]). The 
dashed and dotted curves show the observational fits advocated by Gebhardt et al. (2000) and 
Farrarese & Merritt (2000), respectively. 

Figure 2. The ratio /ig of black hole mass to host mass plotted as a function of the velocity 
dispersion a of the host galaxy. The solid curve shows the prediction of this paper. The data 
points (adapted from Gebhardt et al. 2000) exhibit considerable scatter, but their mean value 
is in reasonable agreement with theoretical expectations; the various symbols represent ellipticals 
(circles), SO galaxies (squares), and spirals (triangles). 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



